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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

SOME FLUTTER EXPERIMENTS AT A MACH NUMBER OF 1.3 ON 
CANTILEVER WINGS WITH TUBULAR AND CLOSED 
BODIES AT THE TIPS 

By John Locke McCarty and W. J. Tuovlla 
SUMMARY 


Flutter tests of 39 cantilever wings with tubular and closed bodies 
on the tips have been made in a small intermittent-flow supersonic wind 
tunnel where the testing technique involved changing the structural param- 
eters so that the wing-body combinations would flutter at the tunnel Mach 
number of 1 . 3 . Some effects of mass moment of inertia of the tube about 
the elastic axis, tube center-of -gravity location, and mass flow through 
the tube on the flutter-speed coefficient were studied. The experimental 
results are compared with a theory using two-dimensional air forces on the 
wing and tube forces derived from a modified slender-body theory. 

The calculated flutter-speed coefficients are in reasonably good 
agreement with the experimental results except at bending to torsion fre- 
quency ratios near 1.0 where the calculated results are much lower than 
the experimental, and at frequency ratios greater than 1.3* where the cal- 
culated results are higher than the experimental. No real flutter-speed 
coefficients could be calculated for the configurations which consisted 
of wings and large closed bodies; however, calculated results were obtained 
for the very slender bodies, on which the body air forces were omitted. 
Divergence calculations were made for two wing-tube combinations and the 
results were unconservative by a wide margin. 


INTRODUCTION 


The aeroelastic phenomena of airplanes and missiles with ram-jet 
propulsion systems located at the wing tips may be influenced to a large 
extent by the flow through the propulsion unit. These effects have 
received little experimental or theoretical attention. Some theoreti- 
cal work, however, has been done on the unsteady forces on slender bodies 
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which would be required in the analyses of these phenomena. Refs. 1 
and 2, for instance, treat the problem of unsteady flow about slender 
closed bodies at supersonic speeds by means of linearized potential -flow 
theory. Some experimental results at subsonic speeds on the flutter of 
open and closed bodies mounted on slender struts are presented in refer- 
ence 3 . These results are compared with an analysis of the flutter and 
divergence of such bodies based on slender-body theory. No known experi- 
mental flutter data at supersonic speeds are available for wings with open 
or closed bodies mounted at the tips . 

This paper presents the results of some tests on the flutter of 
cantilever wings at supersonic speed with tubular bodies, which simulate 
"cold" ram jets at the tips; the effect of internal air flow was inves- 
tigated qualitatively by testing models with solid bodies of revolution 
and comparing the results with those obtained with tubular bodies having 
various amounts of internal flow. The tests were made in an intermittent- 
flow supersonic wind tunnel at a Mach number of 1.3. The experimental 
results are compared with those of an analytical development which con- 
siders two-dimensional supersonic wing forces and tube forces derived 
from slender-body theory. 


SYMBOLS 


b 

f 


f 

a 


8h 

g a 

( I a)w 
( X a)t . 


I 


l 


w 


root semichord, ft 

flutter frequency, cps 

first bending frequency, cps 

first uncoupled torsion frequency, cps 

first bending damping coefficient 
first torsion damping coefficient 

mass moment of inertia of wing about elastic axis, slug-ft^ 
per foot span 

mass moment of inertia of tube about elastic axis, slug-ft^ 
per foot diameter 

tube length, ft 

wing length, in. (see fig. 1(a)) 
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t 
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X 1 


z 


(V‘)w 

A 

P 


average mass of tube, slugs per foot diameter 
mass of -wing, slugs per foot span 


sq uar e of nondimensional radius of gyration of tube about 


elastic axis. 



square of nondimensional 
elastic axis, at root. 


radius of gyration of wing about 

C I a)w 

v > 2 


tube internal radius, in. 

elastic axis position in percent tube length, measured 
from leading edge of tube 

wing thickness, in. 

flutter speed, ft/sec 

divergence speed, ft/sec 


elastic -axis position in percent chord, measured from 
leading edge of wing 


wing center-of -gravity position in percent chord, measured 
from leading edge of wing 


tube center-of -gravity position in percent tube length from 
elastic axis, measured positive behind elastic axis 


tube mass density ratio parameter. 


4m+ 


wing mass density ratio parameter at root. 


gw 

«pb 2 


wing taper ratio, Tip chord 

Root chord 

mass density of air in test section, slugs/cu ft 
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V/b ^ 


flutter-speed coefficient 


MODELS AND TEST METHODS 
Models 


The wings of the wing-tube combinations used in these flutter tests 
were made of either magnesium or wood. All magnesium wings were rec- 
tangular in plan form with a chord of 3 inches and a thickness of 
0.064 inch with the exception of one wing which was O.O33 inch thick. 

The wooden wings were constructed of spruce or birch, laminated or non- 
lamlnated, with both rectangular and tapered plan forms; the rectangular 
wings had chords of approximately 2.50 inches, and the tapered wings 
had a taper ratio ranging from O.373 to 0. 502, with a root chord ranging 
from 2.85 to 3.00 inches. All wings were of constant thickness and were 
hexagonal in cross section (see fig. 1(a)). 

Balsa-wood tubes of l/l6-inch wall thickness coated with glass cloth 
and aluminum tubes of 0.005-0. 010 -inch wall thickness were mounted at the 
tips of the various wings. A general description of the wing-tube com- 
binations is given in the sketch in figure l(a). This figure illustrates 
the manner in which the tubes were mounted on the wing. The wings were 
extended through the tubes to provide a rigid connection between the two, 
but the wing tips were considered to be located where the wings entered 
the tubes. Strips of lead were taped to the tube to change the mass and 
inertia of the tube or to give a desired center -of -gravity location. 

The tubes in general were .7 inches in length with an inside diameter 
of 1 inch. Three of the tubes from which flutter data were obtained, 
however, were O.75 inch in diameter and others were less than 7 inches 
in length. 

Diffuser cones of various sizes and shapes, as illustrated in 
figure l(b), were installed in the noses of 8 tubes in order to investi- 
gate the effects of inlet shape and mass flow through the tube on the 
flutter-speed coefficient. Flutter data were also obtained from three 
closed bodies of revolution, 0.55 _ i n ch maximum outside radius, simulating 
external stores on wing tips (see fig. 1(c)) and from two very slender 
bodies of total mass and inertia equivalent to the tubes (see fig. 1(d)). 
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Same of the geometric, inertial, and structural properties of the 
combinations are listed in table 1. Included in table 1 are the first 
bending mode frequencies and the first torsional mode frequencies. The 
first bending mode frequency was taken as the coupled first bending mode 
frequency, which was obtained by exciting the wing at the elastic axis 
and recording the oscillations . The first torsional mode frequency was 
determined by exciting the model in torsion while supporting the wing 
tip at the elastic axis . Structural damping values were obtained from 
oscillograph records of both frequencies . 


Test Methods 

The models were mounted at zero angle of attack as cantilevers on 
an injector and were tested at a Mach number of 1. 3 in an intermittent- 
flow supersonic wind tunnel having a 9- by 18-inch test section. The 
testing technique involved injecting the models into and retracting them 
from the tunnel while the flow was steady at a Mach number of 1.3. This 
technique was used to avoid any possible flutter that might occur during 
the tunnel starting and stopping transients. The testing procedure was 
to clamp the wing very short to be sure the wing-tube combination would 
not flutter and then to increase the span in successive t unne l runs until 
flutter occurred. If flutter had not occurred when the wing length 
reached 4 inches, the structural parameters or the location of the tube 
center of gravity were changed rather than further extending the wing in 
order to avoid possible interference effects from the opposite wall of 
the tunnel. 

Most of the structural -parameter changes were brought about by 
reducing the root thickness in order to change the bending stiffness, 
or by slitting the leading and trailing edges at the root, in order to 
change the torsional stiffness . The center of gravity of the tube was 
shifted by taping lead strips to the tube. Structural -parameter changes 
or changes in the center of gravity were continued until the combination 
fluttered. A more detailed description of the test methods and photo- 
graphs of the apparatus can be found in reference 4. 

The root strains, position of the model in the tunnel, and the 
static pressure in the test section were recorded simultaneously by a 
recording oscillograph. 


METHOD OF ANALYSIS 


In the method of analysis employed in this paper two-dimensional 
supersonic air force coefficients from reference 5 were used for the 
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vings, and the aerody nami c forces given in the appendix of reference 3 
were used for the tubes and closed bodies of revolution. 

In calculating the flutter-speed coefficients of the combinations, 
use was made of the ratio of the experimental first natural bending fre- 
quency to the torsion frequency and of the measured first bending and 
first torsion damping coefficients. For combinations having rectangular 
wings, linear torsion and parabolic bending mode shapes were assumed, and 
for the combinations having tapered wings the calculated natural mode 
shapes were used. The 70 -percent -span station was taken as representative 
of the entire wing on all models except for some tapered wings having high 
ratios of bending to torsion frequency, where it was necessary to integrate 
the wing parameters along the span in order to obtain a flutter solution. 

Figure 2 shows the relationship between the theoretical flutter-speed 
coefficients, V/boc^ obtained from both a spanwise -integration and a 
representative -section analysis and the ratio of first natural bending 
to torsion frequency for the tapered wing model 20 . • The two 

methods of calculation are in good agreement and it may be noted that 
the spanwise -integration method serves to extend the flutter solution 
to a slightly higher ratio of natural bending to torsion frequency, in 
the sense that it predicts flutter for values of this ratio at which the 
other method no longer predicts flutter (at least in the first mode). 

The aerodynamic forces on the tubes and bodies were based upon a 
modified slender-body theory as in reference 3. This modified theory 
disregards the concentrated forces on the tail section of the tube on 
the premise that both the external and internal flows leave the trailing 
edge of the tube tangentially and are not realined with the free stream. 

The calculations of the flutter-speed coefficients for combinations 
having a diffuser cone did not consider the aerodynamics of the cone, 
and were based on the assumption that the cone did not impede the flow 
through the tube. For the calculations of the two combinations having 
very slender bodies, the body air forces were omitted. 

The divergence velocities were calculated for two wing-tube com- 
binations by equating the aerodynamic moment on the tube (for 00 = 0), 
as taken from reference 3> to the restoring moment applied to the tube by 
the wing. Since the elastic axis on all models was at the midchord of 
the wing and since the linearized two-dimensional supersonic theory 
predicts zero moment about the midchord for these wings, only the tubes 
contributed to the moment in these calculations. 
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RESULTS AND DISCUSSION 
Presentation of Data 


All combinations from which either flutter or divergence data were 
obtained are listed in table 1. This table is divided for convenience 
into four sections: flutter of open tubes, flutter of tubes with diffuser 
cones, flutter of closed bodies, and divergence of open tubes. In the 
table are listed the significant properties of the combinations and the 
experimental and analytical results . 


Experimental Results 

As previously mentioned, the tunnel conditions were fixed and it 
was necessary to change the properties of the combinations to obtain 
flutter. Inasmuch as It is generally difficult in experiments to change 
one parameter without changing others, it is difficult to obtain the 
effect of variations in a single parameter on the quantity of interest, 
in this instance, the flutter-speed coefficient. Consequently, in the 
following paragraphs the effects of various parameters are discussed on 
the basis of a comparison of those combinations which varied in essen- 
tially only one parameter. In this way, the effect of wing taper ratio, 
tube inertia, mass flow through the tube, and center-of -gravity location 
on the flutter-speed coefficient could be noted, although some of the 
conclusions had to be based on comparisons of only two combinations . It 
should be noted that in all of these tests the flutter speed and the tube 
length were constants, so that changes in the flutter-speed coefficient 
reflect changes in the torsional frequency which caused the given com- 
bination to flutter. 

Tube inertia .- The results for combinations 13 and 14 indicate that 
the flutter-speed coefficient increases with increased mass moment of 
inertia of the tube about the elastic axis. In other words, although 
increasing the mass moment of inertia of the tube actually would tend to 
decrease the flutter speed as a result of the reduced torsional frequency, 
the decrease in the flutter speed is apparently less than that in the 
torsional frequency, so that the flutter-speed coefficient is increased. 

Mass flow through tube .- Results for combinations 13 and 30, 27 and 7, 
and 26 and 6 indicate that the flutter-speed coefficient increases with 
decreasing mass flow through the tube; eliminating the internal flow 
entirely increases the flutter-speed coefficient further as shown by the 
results for combinations 37 and 22, 36 and 14, and 35 and 19. 
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Center-of -gravity location .- The effect of the center -of -gravity 
location on the flutter-speed, coefficient vas observed directly since 
the flutter condition in some cases vas approached by shifting the 
center of gravity of the tube rearvard. Thus the flutter-speed coeffi- 
cient is increased by moving the center of gravity of the tube forvard. 

Frequency ratio .- No direct experimental relationship between the 
frequency ratio and the flutter-speed coefficient can be obtained from 
a comparison of combinations listed in table 1. The theoretical vari- 
ation of the flutter-speed coefficient with frequency ratio is, however, 
illustrated in figure 2 for one of the combinations tested. This figure 
is typical of combinations presented in this report. 

It may be seen that the flutter-speed coefficient decreases to a 
minimum value at a frequency ratio near 1.0 and increases rapidly there- 
after with an increase in frequency ratio. A comparison of figures 2 
and 3 indicates that the experiments follow a similar trend, but the 
decrease in the experimental flutter-speed coefficient in the frequency 
ratio region of 1.0 is less pronounced than that indicated by the curve 
of calculated flutter-speed coefficients. In fact, the experimental 
flutter-speed coefficient is almost independent of the bending to torsion 
frequency ratio. 

Divergence ♦ - Wing -tube combinations 38 and 59 diverged before a 
flutter condition could be reached because their forward center -of -gravity 
location tended to increase the flutter-speed coefficient, whereas the 
forward tube location and the thick wings tended to decrease the diver- 
gence speed. 


Comparison of Theory With Experiment 

Flutter calculations were made for all combinations which fluttered. 
Solutions were obtained for the wing -tube combinations and also for the 
combinations of wings and the very slender bodies on which the air forces 
were neglected. No solutions could be obtained, however, for the combi- 
nations of wings and the closed bodies on which the body air forces were 
included. These results indicate that the air forces used on the closed 
bodies are stabilizing inasmuch as these wing — closed-body combinations 
appear to be stable theoretically. 

Flutter-speed coefficient .- A comparison of the calculated and 
experimental results on the basis of the flutter-speed coefficient 
V/bcjo a is made in figure 3* In this figure, the ratio of the experi- 
mental flutter-speed coefficient to the theoretical flutter-speed 
coefficient is plotted against the ratio of first natural bending to 
torsion frequency . In the lower range of this ratio the calcu- 
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lated flutter-speed coefficient is within 20 percent of the experimental 
flutter-speed coefficient. When the natural frequency ratio is near 1.0, 
however, the experimental flutter-speed coefficient tends to be much 
higher than the calculated flutter-speed coefficient. In the natural 
frequency ratio region of 1.2 the measured flutter-speed coefficient 
again agrees with the calculated flutter-speed coefficient, but as the 
natural frequency ratio increases beyond 1.3 the measured flutter-speed 
coefficient becomes less than the calculated; in other words, the cal- 
culations become unconservative . 

The deviations may be explained by an examination of the theoretical 
flutter curve as presented in figure 2, which indicates a type of resonant 
condition near the frequency ratio of 1.0 and an asymptotic behavior at 
ratios greater than 1.2. (For other combinations the asymptotic tendency 
occurs at other frequency ratios . ) As stated previously, the experimental 
flutter-speed coefficient is almost independent of the ratio of bending to 
torsion frequency; the experimental flutter-speed coefficients tend to be 
lower when the ratio (Xt^j is near 1.0 than elsewhere, and they tend to 

increase as this ratio increases beyond 1.0, but these variations are 
much less pronounced than those of the calculated flutter-speed coeffi- 
cients. Therefore, the variation of the ratio of experimental to calcu- 
lated flutter-speed coefficient in figure 3 is primarily due to the 
variation in the calculated flutter-speed coefficient. In general, con- 
sidering slight experimental errors and the fact that wing aspect ratio, 
tube inlet shape, tube mass flow, and interference effects were neglected 
in the theoretical analysis, the agreement between the calculations and 
the experiment for tubes is satisfactory except for the regions where 
equals 1.0 and where 0 ^/ 0 ^ is greater than 1.3. 

Flutter frequency .- Figure 4 is a plot of the ratio of the experi- 
mentally obtained ratio oyco a to the calculated ratio a>/cu a against 

the ratio of the natural bending to torsion frequencies. This figure 
indicates that the theory predicts the flutter frequency reasonably well. 

Divergence . - Divergence calculations were made by using the formula 
in reference 3 and the results are very unconservative; the actual diver- 
gence speed being 1430 feet per second, whereas the calculated divergence 
speeds are 2140 and 2470 feet per second for combinations 38 and 39, 
respectively. Apparently, the divergence calculations greatly underes- 
timate the aerodynamic moment. This discrepancy may be due to two reasons. 
The assumed flow conditions at the tail section may be unrealistic and 
large forces may exist there. Also, the chordwise center of pressure of 
the steady aerodynamic forces on the wing, as given by linearized theory, 
is at the midchord but is known to be actually ahead of the midchord 
position, giving an additional moment not included in the analysis. 
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CONCLUDING REMARKS 


Flutter tests of 39 cantilever wings with tubular and closed bodies 
on the tips were made in a small intermittent -flow supersonic wind tunnel 
where the testing technique involved changing the structural parameters 
so that the wing-body combinations would flutter at the tunnel Mach 
number of 1.3* The following conclusions were made: 

1. The experimental results indicate that within the range of model 
parameters the flutter-speed coefficient is increased (other parameters 
remaining fixed) by increasing the mass moment of inertia of the tube 
about the elastic axis, shifting the center of gravity of the tube for- 
ward, and decreasing the mass flow of air through the tube. 

2. The calculated flutter-speed coefficients are in reasonably good 
agreement with the experimental results except at ratios of bending to 
torsion frequency near 1.0 where the results of the calculations tend to 
be conservative by a wide margin and at ratios of bending to torsion fre- 
quency greater than about 1.3 and up to the highest natural frequency 
ratio of 1.43 where the results of the calculations tend to be somewhat 
unconservative . 

3. A plot of the calculated flutter-speed coefficient against the 
ratio of bending to torsion frequency indicates a minimum flutter -speed 
coefficient near the frequency ratio of 1.0 and thereafter a very rapid 
increase with increased frequency ratio. It appears, therefore, that 
the experimental variation of the flutter-speed coefficient with the 
natural frequency ratio is similar to the calculated but not as pronounced. 

4. No real flutter-speed coefficients could be calculated for the 
combinations of wings and closed bodies; however, calculated results 
were obtained for the very slender bodies on which the body air forces 
were omitted. 

3. The calculated flutter frequencies are in good agreement with the 
experimental flutter frequencies. 

0 

6. Two wing-tube combinations diverged before a flutter condition 
could be reached. The results of divergence calculations based on line- 
arized two-dimensional theory were unconservative by a wide margin for 
these two cases . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 6, 1955* 
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(a) Wing -tube combination. 

Figure 1.- Sketches of combinations. 
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Figure 2.- Calculated flutter-speed coefficient plotted against ratio of 

natural bending to torsion frequency for model 20; both the representative- 
section and spanwise -integration analyses are illustrated. 
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Figure 3.- Ratio of experimental to calculated flutter-speed coefficient 
plotted against ratio of natural bending to torsion frequency. 
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Figure 4 .- Ratio of experimental to calculated cd/o^ plotted against 
ratio of natural bending to torsion frequency. 








